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ABSTRACT 

A method is developed for the optimal design of composite links based on dynamic 
performance criteria directly related to structural modal damping and dynamic stiffness. 
An integrated mechanics theory correlates structural composite damping to the param- 
eters of basic composite material systems, laminate parameters, link shape, and modal 
deformations. The inclusion of modal properties allows the selective minimization of vi- 
brations associated with specific modes. Ply angles and fiber volumes are tailored to obtain 
optimal combinations of damping and stiffness. Applications to simple composite links in- 
dicate wide margins for trade-offs and illustrate the importance of various design variables 
to the optimal design. 


INTRODUCTION 

Fiber reinforced composites show great potential as candidate materials for high-speed 
mechanisms and high-speed robotic manipulators. Continuous fiber composites provide 
superior specific modulus and specific strength allowing the fabrication of stiff, strong, and 
light-weight links. In this iiw.' ner, positioning errors associated with the static or dynamic 
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response are readily decreased, since the elastic deflections are reduced, the strength is 
increased, and the mass of the links is also decreased resulting in lower inertial loading. 
Composites also provide the unique capability to tailor the anisotropy of the material to 
meet specific design requirements. For the case of robotic manipulators, tailoring of the 
material may result in robotic arms with superior and more uniform properties within 
their working envelope [1]. 

In addition to high specific modulus and strength, composite materials will incorpo- 
rate higher material damping into the structure, since most polymer- matrix composites 
exhibit higher damping factors by an order of magnitude compared to commonly used 
metallic alloys. Properly designed composite links can provide higher passive damping, 
and they may further improve accuracy, strength, and fatigue endurance by attenuating 
undesirable elasto-dynamic phenomena such as overshooting, long settling times and struc- 
tural resonances. Recent work [2,3] has shown that structural composite damping is highly 
tailorable, and is primarily controlled by constituent parameters (fiber/matrix properties, 
volume fraction), laminate parameters (ply angles/thicknesses, stacking sequence), and the 
deformation of the articulated structure (mode shapes). The same study has also demon- 
strated that any increase in damping usually results in decreased stiffness and possibly 
strength, therefore, any tailoring of the composite material for optimal damped response 
will be based on trade-offs between damping, stiffness, and strength. In order to realize 
full benefits from the damping capacity of composite materials, integrated methodologies 
for the optimal design of composite links should be developed entailing: (1) capability for 
tailoring the basic composite materials and/or laminate for maximum damping, stiffness, 
and strength; (2) dynamic design criteria based, if possible, on the performance of the 
articulated structure; (3) three-dimensional modeling capacity, as the vibrational response 
of robotic manipulators and mechanisms is in general three-dimensional. The current pa- 
per presents the development of such a formal design method allowing the tailoring of 
composite links for optimal elasto-dynamic performance. 

Limited research has been reported on the optimal design of composite links includ- 


2 



ing damping as design criterion. The flexural damping capacity of laminated beams [4] 
and tubular composite manipulator links [5,6] has been optimized with bounds on design 
parameters and bending rigidity. Ply angles, ply thicknesses, fiber volume fractions, and 
cross-sectional dimensions were utilized as design parameters. The present optimal design 
method is targeted to three-dimensional kinematic configurations. It implements crite- 
ria based on global structural dynamic properties rather than on local properties, and it 
readily incorporates micromechanics for all properties of the composite material, including 
damping. The method can be used to perform each or combinations of the following func- 
tions: (1) prediction of optimal laminate configurations (ply angles and ply thicknesses); 
(2) optimal tailoring of the composite material in terms of volume fraction, matrix and 
fiber properties; and (3) prediction of optimal link shapes. The method provides the unique 
feature to model and optimize the structural damping of individual vibration modes, there- 
fore, superior tailoring of the composite material may result by selective minimization of 
positioning errors associated with the more predominant modes. In view of the later 
feature, the design method may be readily combined with active vibration control, such 
that the elasto-dynamic effects of higher modes are passively minimized by tailoring the 
composite links, while vibrations from lower modes are actively suppressed. 

The performance criteria used in this work are directly related to the damping, static 
deflections, and dynamic deflections of the composite structure. These properties are 
modeled using finite element analysis. Composite damping simulation capabilities have 
been incorporated into a three-dimensional triangular plate element. The optimal design 
problem is formulated in standard non-linear programming form and is numerically solved 
with the feasible directions method. For a preliminary assessment regarding the merits of 
the method, the optimal design of simple links is investigated. 

NOMENCLATURE 

[M 3x3 Extensional laminate damping matrix. 

[A s ] 3x3 Extensional laminate stiffness matrix. 

[B] Strain shape function matrix. 
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Global and modal damping matrices. 

3x3 Coupling laminate damping matrix. 
3x3 Coupling laminate stiffness matrix. 
3x3 Flexural laminate damping matrix. 
3x3 Flexural laminate stiffness matrix. 
Normal modulus. 

On-axis and off-axis ply stiffness matrices. 
Shear modulus. 

Volume ratio. 

Global and modal stiffness matrices. 

Global and modal mass matrices. 

Number of plies in a laminate. 

Modal vector. 

Time. 

Ply thickness. 

Displacement vector. 

Undamped dynamic amplitude. 

Deflections along x,y,z structural axes. 

Stored and dissipated specific strain energies. 
Stored and dissipated strain energies. 
Distance from laminate mid-surface. 
Engineering strain. 

Fiber orientation angle. 

Curvature. 

Poisson’s ratio. 

Stress. 

Specific damping capacity. 


f Bold characters indicate vectors 
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[*] 

Modal matrix. 


Rotations about x,y,z structural axes, 

a? 

Angular velocity. 

Subscripts 

6 

Ply bottom. 

c 

Ply (off-axis). 

e 

Finite element. 

/ 

Fiber. 

i 

Node. 

l 

Ply (on-axis). 

L 

Laminate. 

m 

Matrix. 

n 

n-th mode. 

rp 

Ply interface. 

t 

Ply top. 

Superscripts 

d 

Dynamic. 

L 

Lower. 

o 

Laminate mid-surface. 

s 

Static. 

U 

Upper. 

Direction 

11 

Normal longitudinal. 

22 

Normal in-plane transverse. 

12 

Shear in-plane. 

XX 

In-plane normal, off-axis, x-direction. 

yy 

In-plane normal, off-axis, y-direction. 
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ss 


In-plane shear, off-axis, y-direction. 


COMPOSITE STRUCTURAL DAMPING 

In this section an integrated theory for modeling structural composite damping is de- 
scribed. Analogous theories were utilized for other mechanical properties [7]. The damping 
of the composite link is simulated based on an upward synthesis scheme. First, the on-axis 
damping of each composite material (damping along the material axes) is synthesized from 
the constituent fiber/matrix properties and volume fractions. Then, the damping capac- 
ities of plies subjected to off-axis loading (off-axis plies) are calculated, and subsequently 
the local laminate damping capacity is formulated. Integration of local laminate damping 
over the whole structure provides the global structural damping capacity at a particular 
deformation state. Finite element discretization is utilized for the numerical integration 
of local laminate damping. The same finite element mesh is also utilized to simulate the 
static and dynamic response. 

On- Axis Composite Damping 

A simplified micromechanics theory [2] has been utilized for the synthesis of hysteretic 
on-axis composite damping. The micromechanics equations for on-axis damping are pre- 
sented in the Appendix. The specific strain energy released per cycle due to in-plane 
on-axis cyclic loading is: 

Swi = ^i} T [i>i}{ei} ( 1 ) 

Off-Axis Composite Damping 

Off-axis composite specific damping capacities (SDC’s) are synthesized from the on- 
axis properties and the ply angle. The transformation law is given in the Appendix. The 
off-axis damping matrix [r/> c ] is fully populated due to in-plane coupling between extension 
and shear. As has been demonstrated in ref. [2] coupling provides an additional physical 
mechanism for strain energy dissipation, and consequently, higher damping. The specific 
strain energy released per cycle during cyclic in-plane off-axis loading is, 

= \{<7c} T [i>c}{zc} ( 2 ) 
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Laminate Damping 

The laminate damping is ensembled from the damping contributions of each off-axis 
ply. As shown in the Appendix and ref. [3], the damping capacit 3 r of a general com- 
posite laminate is formally described by three 3x3 fully-populated matrices: (1) the 
extensional damping matrix [Ad] associated with extensional deformations only; (2) the 
flexural damping matrix [Dd] associated with bending and torsional deformations; and 
(3) the extension-flexure coupling damping matrix [Cd] associated with coupled extension- 
flexure deformations. The coupling damping matrix [Cd] is identically zero for symmetric 
laminates. In addition to hysteretic ply damping, the laminate damping includes also 
contributions from hysteretic interlaminar strain energy dissipated in the interply matrix 
layer, as a result of the interlaminar shear stress developed between two adjacent plies 
in order to restrict their relative local movement (scissors effect). The dissipated strain 
energy through-the- thickness of the laminate per unit area per cycle is: 

Wi = i{ £ ° *} 

As seen in eq. (3), the laminate damping is a local structural property, depending on the 
local deformation, the basic constituent composite materials, and laminate configuration. 
Correlation with experimental measurements has yielded excellent agreement [3]. 

Structural Modal Damping 

The local laminate damping provides rough estimates regarding the global damping 
capacity of simple structural components and loading conditions where the local strain and 
curvature remain mostly uniform over the structure. However, most composite structures 
of practical interest would exhibit more complicated and coupled deformations, and their 
global damping capacity would be the integrated contribution of local damping over the 
whole volume at a particular deformation shape. 

Since the elasto-dynamic response of any structure is a linear combination of individual 
vibration modes, structural damping associated with each mode would provide a natural 
and meaningful measure regarding the overall damping capacity of the composite structure. 


[Ad] [Cd] 
[Cd] [D d ] 


e 

K 


( 3 ) 
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In the remaining sections, the global damping capacity of the structure for a particular 
mode shape will be identified as structural modal damping. In this manner, an array 
of modal damping factors is generated, each one representing the damping capacity of 
the structure vibrating solely at the respective mode shape. The modal specific damping 
capacity (SDC) associated to the n-th vibration mode would be: 




J v Sw n dV 
Jy W n d\ 


(4) 


where Sw n and w n represent the dissipated and stored specific strain energies per cycle. 
Utilizing eq. (3), the modal damping of plate/shell composite structures becomes, 
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The integrations in eq. (5) are performed numerically using finite element discretiza- 
tion. Such a finite element scheme has been developed for an isoparametric triangular 
plate element with 3 nodes, 6 degrees of freedom per node Ui = {u x , u yt) u z > <f> x , <j> yti <^ 2 }, 
linear shape functions for membrane deflections, and cubic shape functions for flexure. 

Integration is first performed over the area of the triangular finite element. Given the 
shape functions for membrane and flexural deflections, the mid-plane strain and curvature 
at a point is related to the nodal displacement vector, 



= [B}u e 


(6) 


Based on this relation, the damping and strain energies per cycle, over the area of the 
triangular element would be respectively: 


SW e , n 
W e , n 



U e , n T [B} T 
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Both integrals are calculated numerically using three-point quadrature. Once the dissi- 
pated and stored strain energies are calculated for each element, then their summation 
provides the total dissipated and stored energies. The modal SDC of the n-th node 't/’n is 
then assembled according to equation (5). 

^ ( 8 ) 

j= 1 

ELASTO-DYNAMIC RESPONSE 

The dynamic elastic response of composite links is simulated with finite element anal- 
ysis. The plate triangular element mentioned previously is used for this purpose. This type 
of element was preferred in this initial study, because it provides more accurate modeling 
and insight, than the beam finite elements used in previous studies [1,4-6]. In this manner, 
a reliable data-base of results is assembled for future assessments of the accuracy of less 
expensive beam elements. Using finite element discretization, the elasto-dynamic response 
of the composite link may be simulated by the following system of n discrete dynamic 
equations, 

[M\{u] + [C]{u} + [A']{u> = {F(t)} (9) 

The damping matrix is formulated in the modal space. A typical exciting force for 
composite robotic manipulators would involve: surface tractions F a (t), body forces 
inertial forces due to local, centrifugal and Coriolis acceleration Fi(r), and hygrothermal 
forces Fj/Qn(t). 

W)} = {F.(t)} + {*»(<)} + W?)} + {Fh t(()} (10) 

This study neglects the effects of body and hygrothermal forces, hence, is limited to links 
operating in room temperature and dry environment. Hygrothermal effects may be sig- 
nificant in the design of composite links, and the method outlined here has provisions to 
model them, however, the design problem is complicated beyond the scope of the current 
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paper, as the damping and stiffness of the composite depend also on temperature and 
moisture variations [3]. 

The dynamic system in equation (9) is transformed to the m x m modal space through 
the linear transformation [8]: 

M = W{«} (ii) 

where q is the modal displacement vector, and [$] contains the first m normalized un- 
damped eigenvectors. For proportional damping, the transformation yields the reduced 
uncoupled dynamic system, 


m K<?} + kK$} + = {/(<)} 


( 12 ) 


where, 

[to] = [$] T [Mj[$] 

H = ffl T [C]ffl 
M = !*] T M*1 
{/(<)} = W T mt)} 


The decoupled damping matrix is formulated from the first m modal SDC’s from 
equations (5-8). For general load histories, the differential equations are numerically solved. 

Impulse Response 

The response of the composite link to an impact force {F}h(t), where h(t) is the unit 
impulse, would be [8]: 

_ g — Cn^nt 

Uj(t) = > $jn/n «nK»l) (13) 

, TO n U >d,n 

n—1 

where u >d,n is the damped natural angular velocity of the n-th mode, j any nodal displace- 
ment of the dynamic system (9), and {f}h(t) = [$] T {F}h(t). Eq. (13) is rewritten in the 
form: 

u j( t ) = U ln — r == y f sin (ud,nt) (14) 

n= 1 V 1 U 
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Eq. (14) illustrates the effect of modal undamped amplitudes Uj n and modal damping 
ratios on the impulse response. The undamped amplitudes of each mode, 

TTd _ *&jnfn 
^ j.n 

J m n U) n 

are used as performance measures for the dynamic stiffness of composite links in the 
following section. 



OPTIMAL DESIGN METHOD 

Most common transient vibration problems in high-speed manipulation are usually in- 
duced by impulsive excitation. During minimal time (bang-bang) control, the actuators are 
mostly operating at their maximum torque capacities. The minimal time movement of the 
end-effector between two points typically involves an acceleration phase under maximum 
positive torque, immediately followed by a deceleration phase under maximum negative 
torque. More complicated maneuvers of articulated robotic arms with multiple degrees 
of freedom exhibit the same fundamental pattern of switchings between acceleration and 
deceleration. Impact excitations are applied during the initiation of the motion, each 
switching from acceleration to deceleration and at the end of the rigid body motion. 

It is reasonable, therefore, to tailor the composite links for optimal elasto-dynamic 
response in impulsive excitation. As seen in eq. (14), a feasible way to control the elasto- 
dynamic positioning error would be by restraining each undamped modal amplitude within 
acceptable limits. One advantage of using dynamic amplitudes, in addition to static stiff- 
ness, is that the modal mass is also introduced into the design. In this manner, the 
elasto-dynamic vibration will remain within acceptable bounds, while the maximization of 
modal damping will ensure minimal settling times. The static stiffness is also bounded, 
such that, acceptable static deflections under the maximum static loads are ensured. 

Two different optimal design cases are formulated, each one representing different 
design considerations and requirements. Both cases are presented in standard non-linear 
programming form. 

Case 1: 


11 



The present design case aims to maximize the damping capacity of the first mode, in 
conjunction with other constraints. The first mode has the lower frequency and requires 
the higher settling time. Preliminary investigations regarding the vibration of off-axis 
composite beams have shown that the first mode is a bending mode. Case (1) is most 
suitable when active control is infeasible or undesirable, and the only remaining option 
is passive vibration control. In this context, the optimal design problem is formulated as 
follows: 

max(ipi) (16) 


subjected to: 

( i.) Upper and lower bounds on design variables: 

-90.0° <0i< 90.0° 


0.0 < kf ti < 0.70 


(17.1) 

(17.2) 


where 8 l and kf ti are respectively the ply angle and the fiber volume ratio (FVR) of 
the i-th sublaminate. 

( ii.) Constraints on q s static displacements at the tip of the link: 

u* < U*' U j = 1, q$ (18) 

( iii.) Constraints on m undamped modal dynamic amplitudes as provided by eq. (15) 
for qd displacements at the tip of the link: 


U d 

J,n 


< u 


d,U 


J — 1 } •••} Qd) (XTid 


n 


, m 


(19) 


Case 2: 

This case is suitable for manipulators with active vibration controllers. Active con- 
trollers are most effective up to a frequency limit, therefore, they successfully suppress 
the static deflections and the elasto-dynamic deflections from the first vibration modes. 
Passive attenuation of higher modes is then desirable. The present methodology provides 
the capability to design composite links such that the passive damping of higher modes 
is maximized and the respective undamped dynamic amplitudes remain within acceptable 
limits. 
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Assuming that passive vibration control is desirable for modes mi to m, the following 
proposed performance criteria will maximize the modal damping of the mode with the 
higher dynamic amplitude: 


max(ip n ) , m-i < n < m 


where the mode n is selected such that: 


= max{Uj t 


mi ) 


U d ■ 


j = 1 


( 20 ) 


subjected to: 

( i.) Upper and lower bounds on design variables: 


- 90 . 0 ° < 6 t < 90 . 0 ° 

( 21 . 1 ) 

0.0 < k ft i < 0.70 

( 21 . 2 ) 


( ii.) Constraints on undamped modal dynamic amplitudes for q d displacements at the 
tip of the link and for modes m 2 to m: 

U U < Ufjn j = l,...,q d , and n = m J ,...,m (22) 


Design Vector 

In both cases the design vector includes: (1) the ply orientation angle 0; of each 
sublaminate (group of ±0{ plies), and (2) the fiber volume ratios (FVR’s) kf of each 
sublaminate. By considering only composite material parameters, the current study is 
mainly focused on the tailoring of composite links. The effect of shape optimization and 
its interaction with material tailoring will be addressed in a future study. 

Optimization Algorithm 


13 


Both design cases have been formulated in non-linear programming form and are 
numerically solved with the method of feasible directions. The feasible directions algorithm 
is a primal optimization method, and performs a direct search of the feasible design space 
based on first order sensitivity for objective function and constraints. Primal methods 
are more suitable for designs involving computationally expensive constraints, like the 
present one, since they typically require fewer constraint evaluations than the penalty 
transformation methods. The feasible directions method used in this study incorporated 
active constraint set and line-search strategies for improved efficiency. 

CASE STUDIES 

Assumptions 

Application studies were performed regarding the optimal tailoring of a plate compos- 
ite link 6 in long, 1 in wide, and 0.2 in thick. This simple composite link was selected in 
order to obtain an initial assessment of the merits of the present design methodology and 
to identify the important design parameters. A finite element mesh of 55 nodes and 80 
plate triangular elements was used as shown in Figure 1. The link was assumed clamped at 
one end (nodes 1 to 5) and free at the other (nodes 51 to 55). The composite material was 
high-modulus treated graphite (HM-S) fiber in an epoxy matrix. The mechanical proper- 
ties of the matrix and the fiber are shown in Table 1. The material system was selected 
to match the 50% FVR HM-S/DX210 graphite/epoxy composite for which experimental 
damping measurements are provided [9]. The damping properties of the fiber and matrix 
for this material system were back-calculated from the reported on- axis damping values 
of the 50% FVR graphite/epoxy composite (see Table 1) based on the micromechanics 
equations in the Appendix. 

The basic laminate configurations are shown in Figure 2. Both configurations are sym- 
metric and consist of regular angle-ply sublaminates (Ok/ — 0k/0k-..)- Configuration (a) 
involves one sublaminate, and configuration (b) three. Regardless the number of sublami- 
nates, the link consists of 20 plies through the thickness, each O.Olin thick, such that each 
ply may be manufactured from 0.005in thick unidirectional preimpregnated tape. Each 
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sublaminate may have a different FVR, consequently, the i-th sublaminate will introduce 
2 design variables into the global design vector: the ply angle and the FVR kfj. 

The static load is a uniform transverse (z-axis) force of 5.0 lb/in applied at the tip of 
the link (nodes 51-55), and the dynamic load a uniform transverse (z-axis) impulse of 0.1 
lb/in applied at the same location. The reference design was taken to be a unidirectional 
composite link of 50% HM-S/epoxy, since a unidirectional composite provides very high 
rigidity but low damping in bending. Unless otherwise stated, the static (z-axis) deflections 
along the tip were restricted to be less that 0.05 in (2.5 times the maximum static deflection 
of the unidirectional link), and each individual modal amplitude of the transient (z-axis) 
response at the middle of the tip (node 53) to be less than 1.87 in (2.0 times the maximum 
amplitude for the unidirectional link). 

Case 1 

Three different subcase studies were performed for this type of problem for: (a) one 
sublaminate with fixed 50 % FVR, (b) three sublaminates of 50 % FVR, and (c) three 
sublaminates with varying FVR’s. 

Subcase la. The current link was assumed to consist of one sublaminate with fixed 
fiber volume ratio. This is the simplest tailoring that can be performed on the composite 
link, and the design vector involves only the ply angle 0 1 . The optimal design is shown 
in Table 2. The optimal, d=27° angle-ply symmetric laminate provides 4.5 times higher 
damping for the fundamental mode with respect to the reference unidirectional link and 
4.2 times higher damping for the second bending mode. The trade-offs for these damping 
improvements are higher static compliance by 2.5 times, higher dynamic amplitudes by 
1.5 times and lower natural frequencies. 

Subcase lb. The current subcase involves three sublaminates of the same basic 
material system. Sublaminates 1 and 2 have 4 plies, and sublaminate 3 has 2 plies. The 
design vector now involves three ply angles. The optimal design is shown in Table 2. The 
utilization of three sublaminates has not improved the optimal damping capacity of the 
fundamental mode compared to the optimal design of Subcase la. The resultant optimal 
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ply angles demonstrate that the outer sublaminates are critical to the damping and stiffness 
of the composite link, while the inner sublaminate had a rather insignificant effect. Such 
type of arrangement was indeed expected, because the first bending mode controls the 
design. However, the ±5° plies in sublaminate 3 provide higher in-plane stiffness in x- 
direction than the optimal design of Subcase la. The effect of higher in-plane stiffness is 
indicated by the higher value of the second natural frequency. 

Subcase lc. In this subcase, the fiber volume ratio of each sublaminate system is 
incorporated into the design vector, which now consists of 3 ply angles and 3 fiber volume 
ratios. As seen in Table 2, the introduction of FVR's as design variables has produced a 
sandwich type of laminate, in that, the fiber volume ratio of the outer sublaminate has 
reached its upper bound. The outer sublaminate provides more bending stiffness but less 
damping. The FVR's of the other two sublaminates have been slightly changed. This 
laminate arrangement has produced almost 8% higher SDC for the first modal damping 
and higher SDC for the second bending mode, without any additional penalties on the 
static stiffness, dynamic amplitudes, and natural frequencies. 

All subcases demonstrate wide margins for trade-offs between various modal damping 
capacities, static and dynamic stiffness, and natural frequencies. Additional insight into 
the dynamics of the problem is provided by Figures 3 and 4, which illustrate the mode 
shapes for the unidirectional reference design and the optimum design of Subcase lc . The 
bending-torsion coupling is minimal. The optimum design has also switched the order of 
modes and has affected the SDC’s of the torsional modes. Torsional modes may also induce 
elasto-dynamic deflections at the end-effector of robotic manipulators. Although the cur- 
rent applications are restricted to transverse bending deflections, the torsional deflections 
may be significant for manipulators with multiple links. 

Case 2 

A link with the previous laminate configuration of three sublaminates and three differ- 
ent material systems was optimized. The basic assumption of the current case is that the 
fundamental mode will be eliminated with active vibration control, hence, the optimum 
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tailoring is now solely focused on modes 2 to 5. In this context, the damping capacity 
of the most predominant mode between modes 2 and 5 was maximized. The imposed 
upper bound on the modal dynamic amplitudes of modes 2 to 5 was 0.29 in (2 times the 
maximum modal amplitude for modes 2 to 5 of the unidirectional link). The optimum 
design is shown in Table 3. The mode with the higher undamped amplitude is mode 3 
(2nd bending), while at the reference design the most predominant mode was mode 4 
(2nd bending). The respective modal damping has increased by a factor of 5, while the 
respective modal amplitude has been increased by a factor of 2. The damping capacity of 
the fundamental mode 1 (1st bending) has been also drastically increased, although was 
not incorporated into the design criteria. The current case illustrates that the utilization 
of active controllers in composite robotic manipulators, will not only eliminate the elasto- 
dynamic effects of the lower modes, but it will also facilitate the design of composite links 
with better passive vibration suppression for the higher modes. 

The transient responses of the unidirectional and optimum (subcase lc) designs are 
shown in Figure 5. Each graph presents the tip vibration (node 53) of the respective 
design induced by the z-axis impulsive load of 0.1 lb/in applied at the tip of the link. The 
superiority of the optimum design is clearly illustrated. Initially, the dynamic deflections of 
the optimum design were higher than the dynamic deflections of the unidirectional design, 
however, because of the higher damping the deflections of the optimum design were quickly 
reduced to lower values compared to the unidirectional design. At 0.6 secs, the deflections 
of the optimum design were virtually reduced to 40 % of the respective deflections of the 
unidirectional link. 


SUMMARY 

A formal design method was presented which automates the design of composite 
links for optimal elasto-dynamic response. The method was based on static and dynamic 
performance criteria. The proposed dynamic design criteria were directly related to the 
damping and amplitudes of selected vibration modes. The elasto-dynamic response was 
simulated with finite element analysis. Simplified, design oriented micromechanics and 
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laminate mechanics models for composite damping correlated the performance criteria 
with laminate and composite material parameters. The present method can be utilized for 
basic material, laminate, and structural shape optimization. Two distinct design problems 
were formulated aiming to manipulators with and without active controllers. 

The applications illustrate the effect of ply angles, multiple angle-ply sublaminates, 
and fiber volume fractions on the optimal design. Based on the case studies, the following 
conclusions are summarized: 

1. The resultant optimal laminate configurations and composite material systems 
significantly improved the damping of bending modes by factors of 5 with respect to the 
reference unidirectional design, while the undamped dynamic amplitudes only degraded 
by factors of 1.5. Wide margins for trade-offs exist between damping, and static/dynamic 
stiffness. These trade-offs can be achieved without changing the geometric shape of the 
links, and consequently, their mass. 

2. Ply angles were the most critical design parameters. 

3. Laminates with single or multiple angle-ply sublaminates produced optimal designs 
of comparable performance. 

4. The optimization of fiber volume fractions has produced small improvements in 
the damping of bending modes. 

5. Dynamic properties associated with torsional modes where also found sensitive to 
variations of laminate and composite material parameters. The actual effect of torsional 
modes on the end-effector deflections is yet unclear. Minimization of the torsional dy- 
namic effects will be better accomplished with design criteria based on the elasto-dynamic 
response of the articulated robotic arm, rather than on individual links. 
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APPENDIX 


The present Appendix summarizes the micromechanics and laminate mechanics theo- 
ries for composite damping implemented in this work [2,3]. Analogous micromechanics and 
laminate mechanics theories [7] have been utilized for the remaining mechanical properties 
of polymer matrix composites. 

On-axis Composite damping: 


Longitudinal normal ply SDC tpm: 


I / i E | I i •L'm 

Vlll — V’ /II "-/~ d" Vmn^/7: 

■U'm &i\\ 


Transverse normal ply SDC’s 4>i22'- 


4>122 = 4'f22\fk~S^~ + ^mn(l “ \/kf)^P- 
&f 22 Zm 


where, 


E 22 — (1 — \/ r kf)E rn + 


VKfEr 


1 ^7) 


In-plane shear ply SDC ^ 112 ' 


12 


i'll 2 = i’mV^f-TT*- + ^m.(l - y/kf)~L 

07/12 Or m 


where, 


£12 — (1 — x/kf)G m + 


■JkjG, 


1 - V*7(i- <fc) 


Off-Axis Composite Damping: 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 


where, 
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(A 8) 
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Laminate Damping: 

The SDC is expressed in the form of extensional, coupling, and flexural 3x3 damping 
matrices. 


Nl 


Nt —1 


[Ad] = X (zt - + X! H ^ S ']*#» 
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[£*] = X ~ Z l)i[Ec]i[i>c)i -T X ZrpMSU^mji (All) 

i=l i=l 
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1=1 1=1 

[i'm] is the matrix damping, Hj the interply distortion energy coefficient [7], and [E c ] the 
off-axis ply stiffness matrix, 
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Table 1 Mechanical properties of the HM- S/epoxy system. 


Epoxy 

HM-S Graphite 

50% HM-S/Epoxy 

E m = 0.500 msi 

Efu =- 55.0 msi 

= 0.45 % 

G m = 0.185 msi 

Ef22 = 0.9 msi 

^122 = 4.22 % 

i' mn = 10.30 % 

G f\2 — 1-1 nisi 

■0112 = 7.05 % 

i'm, =11-75 % 

1^/12 = 0.20 



V’/n = 0.4 % 



V 722 = 0.4 % 



ifrfi 2 = 0.4 % 




Table 2 Optimum designs for Case 1. 


Property 

Unidirectional 

Optimum 
(Case la) 

Optimum 
(Case lb) 

Optimum 
(Case lc) 

Ply Angles, degrees 

9i 

0.0 

27.00 

27.16 

30.24 

@2 

0.0 

- 

27.30 

30.49 

0 3 

0.0 

- 

5.21 

29.76 

Fiber volume ratios 

kf i 

0.50 

- 

- 

0.69 

k f2 

0.50 

- 

- 

0.53 

fc /3 

0.50 

- 

- 

0.50 

Modal SDC’s, (%) 

Pi 

0.571 

2.60 

2.5S 

2.78 

'll' 2 

5.566 

3.90 

2.41 

3. 87 

Pz 

3.426 

2.74 

2.71 

2.81 

P\ 

0.656 

0.82 

0.82 

0.63 

p5 

4.540 

2.22 

2.20 

2.18 

Max. static deflection (tip), (in) 

0.0199 

0.0500 

0.0500 

0.0499 

Modal amplitudes (tip), (in) 

U d zl 

0.934 

1.536 

1.534 

1.518 

U*2 

0.000 

0.000 

0.000 

0.000 

ui 3 

0.000 

0.244 

0.243 

0.240 

v * 4 

0.145 

0.000 

0.000 

0.000 

ut 5 

0.000 

0.091 

0.091 

0.087 

Natural Frequencies, (Hz) 

h 

383. 

252. 

252. 

248. 

f 2 

1089. 

1360. 

1534. 

1338. 

/ 3 

1565. 

1572. 

1572. 

1558. 

A 

2395. 

2964. 

2968. 

3518. 

h 

3854. 

5029. 

5027. 

5053. 
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Table 3 Optimum design for Case 2. 


Property 

Unidirectional 

Optimum 

Ply Angles, degrees 


0.00 

31.24 

0 2 

0.00 

31.74 

0 3 

0.00 

30.02 

Fiber volume ratios 

fc/i 

0.50 

0.53 

k/ 2 

0.50 

0.52 

h* 

0.50 

0.52 

Modal SDC’s, (%) 

Vi 

0.571 

3.35 

V?2 

5.566 

4.10 

i'z 

3.426 

3.34 


0.656 

0.70 

Vs 

4.540 

2.59 

Max. Static deflection (tip), (in) 

0.0199 

0.0660 

Modal amplitudes (tip), (in) 

VU 

0.934 

1.775 

ui 

0.000 

0.000 


0.000 

0.280 

uU . 

0.145 

0.000 

v d :i 

0.000 

0.101 

Natural Frequencies, (Hz) 

h 

383. 

220. 

} 2 

1089. 

1275. 

h 

1565. 

1383. 

u 

2395. 

3274. 

fs 

3854. 

4488. 
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(B) FINITE ELEMENT MESH. 
FIGURE 1. - TYPICAL COMPOSITE LINK. 


-<V e 1>5' k f1 


(A) ONE SUBLAMINATE. 

(e 1 /-e 1 ) 2 , k (1 

( 02 ^~^ 2 ^ 2 ' k t2 

(0j/-e 3 ), k (3 

(B) THREE SUBLAMINATES. 

FIGURE 2. - TYPICAL LAMINATE CONFIGURATIONS. 
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TIME, sec 

(B) OPTIMALLY TAILORED COMPOSITE LINK, 
CASE (1c). 

FIGURE 5. - TRANSIENT IMPULSE RESPONSE. 
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FIGURE 4. - MODE SHAPES OF THE OPTIMUM DESIGN, CASE (1C) 
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